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ABSTRACT 


Most previous studies identifying Panopea generosa and P. globosa have used non-rigorous 
visual methods as well as older shell measurement techniques. Newer mathematical methods 
based on shell shape variation allow for more accurate identification of clam species, as well 
as modeling of phenotypic differences due to environmental effects in populations in different 
sites. Interspecific shell morphology for two Mexican geoduck clam species was analyzed 
from a total of five sites off both coasts of the Baja California peninsula. In addition, intraspe- 
cific analyses of shell morphology were conducted for one of the species, P. globosa, at four 
sites along its reported distribution. Two approaches were employed for the analyses: a novel 
approach based on radiating lines to characterize shell outlines, and a more traditional ap- 
proach using internal shell landmarks. In general, the novel approach afforded greater fidelity 
in distinguishing inter- and intraspecific variation. Our results from both methods agree with 
original species descriptions, and showed that Bahia Magdalena geoducks are P. globosa, 
thus revealing a wider distribution than previous reports for this species. The outline and 
internal scars were highly discriminant between the two species. Shell shape of P. generosa 
was also less variable than that of P. globosa. Intraspecific analyses of P. globosa shell shape 
suggest an adaptive or phenotypic response to environmental conditions at each site. Our 
results may also be indicative of reproductive isolation between Pacific P. globosa at Bahia 


Magdalena and conspecifics in the Gulf of California. 
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INTRODUCTION 


Mollusks of the genus Panopea include among 
the largest of all deep-burrowing bivalves, with 
congeners occurring worldwide in intertidal 
and subtidal regions in the Pacific and Atlantic 
Oceans, the Mediterranean Sea, as well as off 
the coasts of Australia and New Zealand (Yonge, 
1971). Two Panopea species, P. generosa 
(Gould, 1850) and P. globosa (Dall, 1898), com- 
monly called the “geoduck clam”, occur in coastal 
waters of the Baja California Peninsula in both 
the Pacific Ocean and the Gulf of California. 

Some authors report P. generosa distributions 
from south of Alaska to the Baja California pen- 


insula (Weymouth, 1921; Fitch, 1952; Moore, 
1968; Morris et al., 1980; Goodwin & Peace, 
1989; Coan et al., 2000), whereas others in- 
clude Baja California Sur (SAGARPA, 2007) 
and Gulf of California (DFO, 2000; Moore, 
2001; Feldman et al., 2004; Demere & Scott, 
2006) within its distribution range. In contrast, 
the distribution of P globosa has been reported 
as limited to the northern Gulf of California 
(Keen, 1971; Coan et al., 2000; Aragon-Norie- 
ga et al., 2007; Rocha-Olivares et al., 2010), 
although specimens were recently reported 
from the central eastern side of the Gulf of 
California off the mainland coast of Mexico 
near Guaymas, Sonora (Aragon-Noriega et 
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FIG. 1. Sampling sites of two Panopea species in northwest Mexico. 


al., 2007), as well as the Pacific side of the Baja 
California peninsula in coastal waters of Bahía 
Magdalena (Fig. 1; González-Peláez, Centro de 
Investigaciones Biológicas del Noroeste, México 
2010, pers. comm.). 

Misidentifcation of Panopea species, on the 
other hand, remains a problem since morpho- 
metric plasticity is known to occur in the genus. 
For example, one case that is well documented is 
that of Panopea generosa, which had been mis- 
takenly synonymized with the extinct P. abrupta 
for almost 25 years (Vadopalas et al., 2010). 

Moreover, identification can be difficult be- 
cause bivalve shell morphology is influenced 


by such environmental factors as temperature, 
tidal excursion, wave exposure, water currents 
and sediment type (Costa et al., 2008a). The 
marine environment along the peninsula and 
the Gulf of California is highly variable, with 
seasonal circulation patterns (Zaytsev et al., 
2003; Kessler, 2006; Lavín & Marinone, 2003) 
and differences in bathymetry (Bray, 1988; 
Lluch-Cota et al., 2007) that may influence 
regional morphological changes in metapopu- 
lations (Helenes & Carreño, 1999; Bizzarro, 
2008; Álvarez et al., 2009). 

The relationship between local environmental 
conditions and shell morphology in Panopea 
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TABLE 1. Geoduck (Panopea spp.) sampling sites and number 
of specimens used for interspecific and intraspecific geometric 


morphometric analysis. 


Locality Abbreviation 
San Quintin SQ 
San Felipe SF 
Puerto Peñasco RP 
Guaymas GU 


Bahía Magdalena BM 


are unknown. Given the propensity for vari- 
ability in valve shape due to regional environ- 
mental variation (Fyfe, 1984), characterizing 
the variation in shell morphology between and 
within Panopea species could allow for more 
accurate estimates of population distributions. 
Shell morphology and structural traits are the 
principal characters used in traditional taxon- 
omy of bivalves, including Panopea generosa 
and P. globosa. The shell of adult P. generosa 
is heavy, ponderous, subquadrate, subequilat- 
eral and the valves gape widely (Coan et al., 
2000; Vadopalas et al., 2010). By comparison, 
the shell of adult P globosa is more inflated, 
expanded and rounded in the pedal region, 
with a narrower posterior end and non-parallel 
opposite margins (Dall, 1898). Fyfe (1984) 
described the shell morphology of P. generosa 
as highly variable and taxonomically difficult 
to identify using traditional morphometric data 
(length, depth, and width). Thus, morphological 
differences between congeners of Panopea 
and its conspecifics are often difficult to identify 
by inexperienced observers. 

Because linear measurements can be fairly 
ambiguous (Zelditch et al., 2004), researchers 
have explored geometric methods as an alter- 
native for quantifying and analyzing morpho- 
logical variation. Some of the more promising 
methods include the use of points that fall at 
defined intervals along a curve between two 
distinct landmarks; Marko & Jackson (2001) 
termed these points semi-landmarks. 

Geometric morphometrics in other taxa 
have been shown to be powerful tools for 
solving complex species-level identification 
problems (Matias et al., 2001; Gumiel et al., 
2003; Shipunov & Bateman, 2005). Relatively 
few studies have applied these techniques to 
mollusks. Some focused on general morphol- 
ogy (Ubukata, 2003; Anderson & Roopnarine, 
2005; Carvajal et al., 2005, 2006; Roopnarine 


Location Specimens 
30°23’N, 115°57’°W 30 
31°01'N, 114°49°W 18 
31°18 N, 11333 W 30 
27T 50N, 10°51 30 
24°39’N, 112°03’W 30 


et al., 2008), while others used geometric 
morphometrics to compared external shapes 
among conspecific populations (Kwon et 
al.,1999; Costa et al., 2008b). or phenotypic 
stocks (Marquez et al., 2010). 

To help resolve the difficulty in species 
identification, which is affected by shell shape 
variation, we used two methods of geometric 
morphometric analyses to characterize shell 
shape differences in Panopea spp. from five 
sites on the Baja California peninsula. The first 
method is Procrustes superimposition which, 
to our knowledge, has not yet been used with 
bivalves. This newer method involves generat- 
ing radiating lines, which fan out from a centroid 
and is referred to herein as “Fan-Based.” The 
second method is based on anatomical land- 
marks of the inner shell. 


MATERIALS AND METHODS 
Samples and Sites 


For taxonomic identification, we used the 
original descriptions (Conrad, 1849; Gould, 
1850; Dall, 1898), recent identification keys 
(Keen et al., 1971; Coan et al., 2000) and the 
direct help of Paul Valentich-Scott (Curator of 
Malacology, Santa Barbara Museum of Natural 
History). Using SCUBA, divers collected adult 
P. generosa specimens from San Quintin 
(n = 30) and a total of 108 adult P. globosa 
specimens from San Felipe (n = 18), Puerto 
Peñasco (n = 30), Guaymas (n = 30) and Bahia 
Magdalena (n = 30; Table 1; Fig. 1). 


Imaging and Measurement 
The principles of the FB method are well de- 


scribed elsewhere (e.g., Sheets et al., 2006a; 
Rohlf, 1990; Lawing & Polly, 2009). For this 
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FIG. 2. Digitized landmarks and coordinates located for Panopea generosa (left) and P. globosa (right). 
Where (a) is the fan-based method showing the intercept between the valve contour and the drawing 
aligment fan (48 equal quadrants), and (b) is the internal scar method showing landmarks 1—3 which 
describe nymph shape while landmarks 4—6 describe the deep pallial sinus. 


method, the set of points was defined a priori by 
an arbitrary rule: rays at equal angular intervals 
were used as guides to digitize semi-land- 
marks, with the points to be digitized located 
at the intersection of all radii of a circle and the 
outline curve (Sheets et al., 2006a; Rohlf, 1990; 
Lawing & Polly, 2009). As employed here, the 
FB method involved generating three anatomi- 
cal landmarks and two-dimensional Cartesian 
coordinates of 48 semi-landmarks (Fig. 2a) 
using Makefan6 (Sheets et al., 2006b). We 
used these landmarks to align shapes via Pro- 
crustes superimposition. We used anatomical 
scars for landmarks 1, 2, and 3 (Fig. 2a) (umbo, 
antero-ventral adductor muscle, and postero- 
ventral pallial sinus) to determine centroid size, 
semi-landmarks 4 to 51 to analyze curves and 


shell shape variation (semi-landmarks 4—11 
and 45-51 to define the dorsal region; 12—20 
to define the posterior end; 21—33 to delineate 
the ventral region, and 34—44 to define the 
anterior end). We used marks 52 and 53 as 
a reference scale (tpsDig 2.12, Rohlf, 2008). 
We used Semiland v6 software (Sheets et al., 
2006b) to conduct FB analyses. 

The ILB method analyzed the interior shape 
of each specimen using 11 landmarks (Fig. 
2b). We used landmarks 1-3 to analyze nymph 
shape; 4—6 to analyze pallial sinus depth; 
landmark 7 to characterize the antero-ventral 
adductor scar and landmarks 8—11 to define 
outline curves. As above, we used landmarks 
12 and 13 as a reference scale (Fig. 2b; tpsDig 
2.12, Rohlf, 2008). We translated the interior 
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TABLE 2. Variance explained by the first two principal components between 
Panopea generosa and P. globosa and within P. globosa groups. 


Group 
Outline 
Interspecific P. generosa 
P. globosa 
Intraspecific SF 
PP 
GU 
BM 
Internal landmarks 
Interspecific P. generosa 
P. globosa 
Intraspecific SF 
PP 
GU 
BM 


landmark coordinates to align the centroids 
of each individual. To minimize least-squares 
differences between landmark alignments, we 
used CoordGenéf (Sheets et al., 2006b) to ro- 
tate and scale centroid size to unity, producing 
Procrustes shape coordinates and consensus 
shapes for each group. 


Analyses of Interspecific and Intraspecific Varia- 
tion 


To determine differences in outline (FB 
method) and internal landmarks (ILB method) 
between FP. generosa and P. globosa, we used 
CVAGen6m (Sheets et al., 2006b) to perform 
canonical variates analysis and multivariate 
analysis of variance (CVA-—MANOVA). To as- 
sess within-group shape variation, we used 
PCAGen6n software (Sheets et al., 2006b) to 
perform principal component analysis (PCA). 
To determine which canonical variates (CV) 
were effective discriminators, we employed 
Bartlett's test based on Wilk’s lambda (A), 
which is the within-groups sum of squares 
divided by the total sum of squares. A value of 
A near 0 (zero) indicates high discrimination 
between groups (Zelditch et al., 2004). 

We assessed the statistical significance 
of shape differences using Goodall’s F-test. 
We calculated the full and partial Procrustes 
distances between means (DBM) of two 


PC I (in%) PCIl(in%) Overall (in %) 
52 19 71 
36 30 66 
48 15 53 
35 22 oy 
41 Z5 66 
43 28 71 
36 25 58 
28 17 45 
34 14 48 
oa 19 50 
o 16 51 
oo 21 54 


groups and bootstrapped to estimate the 
variance of this distance using TwoGroup6A 
software (Sheets et al., 2006b). We employed 
Mahalanobis distances (Mahalanobis, 1936) 
to determine the probability that the distance 
between an individual and the mean of the 
group was larger than expected under the null 
model of random variation around the mean of 
each group (Zelditch et al., 2004). 

To visualize the shapes, we used thin-plate 
splines of the first two principal components. 
Thin plate splines provide graphical representa- 
tions of the deformation of a square grid based 
on the difference in position of individual land- 
marks and displacement of the vectors relative 
to a grand mean consensus form (Bookstein, 
1991). 


RESULTS 
Interspecific Differences 


Fan-Based Method: CVA-MANOVA analysis 
of the FB method revealed a single CV with sig- 
nificant differences between species (Bartlett’s 
test, Wilk’s A = 0.0836; p < 0.01; Fig. 3a). Based 
on Mahalanobis distances, all specimens 
(Panopea generosa, n = 30 and P. globosa, n 
= 108) were correctly assigned to species. The 
partial Procrustes DBM between the two spe- 
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FIG. 3. Plot of interspecific canonical analyses of sets of variables and thin-plate 
spline with vectors of outline (a) and internal scars (b), nymph and pallial sinus, 


oval and circle respectively. 
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TABLE 3. Goodall’s F-test and distances between 
means (DBM) of outline and internal landmarks in 
Panopea globosa. BM = Bahia Magdalena; GU = 
Guaymas; SF = San Felipe; PP = Puerto Penasco. 
All comparisons were significant (P < 0.01). 


Comparison F DBM 
Outline 
BM vs GU 16.08 0.0475 
BM vs SF 6.25 0.0344 
BM vs PP 25.45 0.0589 
Guws SF . 2976 0.0440 
Guvs PP 14.94 0.0411 
SF vs PP 8.18 0.0342 
Internal landmarks 
BM vs GU 6.65 0.0540 
BM vs SF 6.35 0.0612 
BM vs PP 11.46 0.0702 
GU vs SF 7.90 0.0636 
GU vs PP 6.96 0.0517 
SF vsiPP 2.69 0.0366 


cies was significant (DBM = 0.065, Goodall’s 
F-test with 900 bootstraps, F = 43.8, p = 0). The 
first two principal components (PC) accounted 
for 71% and 66% of the total variance in P. 
generosa and P. globosa, respectively (Table 
2). Note that the loadings of all other PCs were 
less than 1%; they are not shown in Table 2, as 
they contribute little to the variance. In contrast, 
PC1 and PC2 fan rays contribute to most of the 
variance for the two species. 


Inner Landmark-Based Method: CVA- 
MANOVA analysis of the ILB results revealed a 
single CV with significant differences between 
species (Bartlett’s test, Wilk’s A = 0.0987, p < 
0.01; Fig. 3b). All Panopea generosa and P. 
globosa specimens were correctly assigned to 
species based on Mahalanobis distances. The 
partial Procrustes DBM of internal landmarks 
was significant (DBM = 0.0982, Goodall’s F-test 
with 900 bootstraps, F = 32.72; p = 0). The 
first two PCs accounted for 58 and 45% of the 
total variance in P. generosa and P. globosa, 
respectively (Table 2). 

When we compare the FB and ILB results, the 
thin plate splines of the first two PCs obtained 
by Procrustes superimposition of centroid size 
reveal easily perceivable shape differences in 
both the deformation grids and vector displace- 
ment (Fig. 3a, b). 


TABLE 4. Intraspecific individual assigned with the 
Mahalanobis distance approach among groups 
of Panopea globosa. BM = Bahia Magdalena; 
GU = Guaymas; SF = San Felipe; PP = Puerto 
Penasco. 


BM GUY «im GU 


Outline 
BM 30 0 0 0 
GU 2 28 0 0 
SF 0 0 18 0 
BP 0 1 0 29 
Internal landmarks 
BM 28 1 0 1 
GU 2 27 0 1 
SF 0 1 10 7 
PP 1 2 9 18 


Intraspecific Differences 


Fan-Based Method: Using the FB approach 
we observed three CVs of outline shape 
variation among P. globosa groups with CVA- 
MANOVA (Fig. 4a). Axis 1 (Wilk’s A = 0.0159, 
p < 0.01) separated the Pacific group (Bahia 
Magdalena) from the Gulf of California groups 
(Guaymas, San Felipe, and Puerto Peñasco); 
Axis 2 (Wilk’s A = 0.0965, p < 0.01) separated 
the southern Gulf group (Guaymas) from the 
northern Gulf groups (San Felipe and Puerto 
Peñasco); Axis 3 (Wilk’s A = 0.3620, p < 0.01) 
placed San Felipe and Puerto Peñasco in the 
same group. Mahalanobis distances correctly 
assigned 100% of the specimens to the Bahia 
Magdalena and San Felipe groups. In contrast, 
correct allocation was slightly lower for the 
localities of Guaymas and Puerto Peñasco 
(93 and 97%, respectively). The first two PCs 
account for a high percentage of the total vari- 
ance (58-71%) at all sites (Table 2). Goodall’s 
F-test values and the partial Procrustes DBM 
among P. globosa groups outlines were signifi- 
cant (Table 3). 


Inner Landmark-Based Method: When 
comparing shape variation among internal 
landmarks using the ILB approach in the same 
groups (Fig. 4b), we obtained two CVs: Axis 
1 (Wilk’s A = 0.0800, p < 0.01) discriminated 
the Pacific group (Bahia Magdalena) from the 
Gulf groups (Guaymas, San Felipe, and Puerto 
Peñasco), and Axis 2 (Wilk’s A = 0.3396, p < 
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FIG. 4. Intraspecific canonical analyses of sets of variables plot of shell shape 
variation of Panopea globosa collections in northwest Mexico. Outline (a), 


internal scars (b). 
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0.01) included all Gulf groups. Group assign- 
ment with Mahalanobis distances revealed 
some incorrect assignments of northern Gulf 
geoducks: 37% of Puerto Peñasco specimens 
were grouped incorrectly with San Felipe, 
while 45% of the samples of San Felipe were 
incorrectly grouped with Puerto Peñasco. 
In contrast, > 90% of Bahía Magdalena and 
Guaymas specimens were correctly assigned 
(Table 4). 

With these landmarks, the first two PCs ac- 
counted for 50—54% of the total variance (Table 
2). Goodall’s F-test values and DBM were also 
significant (Table 3). Thin plate splines revealed 
the greatest differences in vector length and 
displacement of landmarks 4—7 and 10 be- 
tween Pacific P. globosa and the three Gulf 
groups (Fig. 4b). 


DISCUSSION 


In this study we used two methods to ana- 
lyze interspecific and intraspecific phenotypic 
variation: a new proposed FB method, for ex- 
ternal contours of the shell, as well as the 
more traditional ILB method for mapping the 
interior of the shells for geometric morphology 
analyses to determine interspecific shell shape 
differences in Mexican geoducks. Our results 
clearly delimited two groups that correspond 
to the two species, Panopea generosa and P. 
globosa, in agreement with the original descrip- 
tions (Gould, 1850; Dall, 1898). The specimens 
from Bahia Magdalena (Pacific coast) grouped 
with P. globosa, corroborating the traditional 
taxonomical procedure. The results indicate a 
much wider distribution for P. globosa than only 
the northern Gulf of California as previously re- 
ported (Coan et al., 2000; Briggs, 1974; Brusca, 
1980; Hastings, 2000; Calderon-Aguilera et al., 
2010; Rocha-Olivares et al., 2010). 


Interspecific Shell Shape Variation 


Prior studies using traditional morphometric 
analyses have shown notable shape differ- 
ences between Panopea generosa and P. glo- 
bosa (Rocha-Olivares et al., 2010). Using three 
linear measurements, these authors found 
that P generosa had a more homogeneous 
shell length compared to P. globosa and they 
successfully discriminated the two species. 
The present study corroborates this finding 
and increases the resolution of inter-specific 
morphological shell variation. 


The ILB method has been used to determine 
shell shape variation in other bivalves. For 
example, Innes & Bates (1999) determined 
contour differences between Mytilus edulis 
and M. trossulus (Wilk’s A = 0.625), whereas 
Roopnarine & Vermeij (2000) reported Wilk’s 
A = 0.25 for internal shape variation between 
Chione cancellata and C. elevata. In contrast, in 
our study Bartlett’s test was highly discriminant 
between the two Panopea species (Wilk’s A 
= 0.08 and 0.09 for FB and ILB, respectively). 

The consensus outline shape of Panopea 
globosa is more inflated, with a narrower pos- 
terior end compared to P. generosa. The umbo 
in P. globosa is near the posterior end, while 
more centrally located in P. generosa. The off- 
center placement of the umbo in P. globosa 
may be related to an anatomical adaptation 
which aids in maintaining valve gape due the 
more inflated shell; Savazzi (1987) discusses 
the hypothesized relationship between shell 
inflation and umbo displacement, and Ubukata 
(2003) observed in pteroids a tendency of shell 
extension in the antero-ventral direction. 

In our study, variation in the inner shell in- 
dicates that Panopea globosa has a shorter 
nymph and smaller ligament (landmarks 1-3), 
but a deeper pallial sinus than P. generosa (land- 
marks 4-6). Stanley (1970, 1979) and Lockwood 
(2004) describe the pallial sinus as the space 
where siphonate bivalves contract the siphon. 
Kondo (1987) suggested that the ratio of pallial 
sinus depth to shell length is closely correlated 
to the maximum burrowing depth. Although we 
found a deeper pallial sinus in P. globosa than 
in P. generosa, we did not measure syphon 
length or burrowing depth thus we are unable 
to corroborate Kondo’s hypothesis. 

The FB and ILB methods showed that the first 
two PCs in Panopea generosa accounted for 
more variability in shell shape than in P. globosa 
(Table 2). To fully corroborate these results in 
a future study, it will be necessary to examine 
specimens of P. generosa from other sites. 


Panopea globosa Intraspecific Shell Shape 
Variation 


Both the FB and ILB methods revealed that 
P. globosa of Bahia Magdalena are morpho- 
logically distinct from the three other localities 
of the Gulf of California. Nevertheless, the FB 
method more clearly distinguished three dif- 
ferent groups, while the more traditional ILB 
method exhibited lower fidelity since it was only 
able to distinguish two groups. 
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Our results suggest a shell shape adapta- 
tion or phenotypic response to environmental 
conditions at each site. The intraspecific CVA 
analyses reveal significant differences among 
Panopea globosa sites, where Pacific P. glo- 
bosa (Bahía Magdalena) have more elongated 
shells but a shorter pallial sinus than those in 
the Gulf of California. Moreover, within the Gulf, 
our analyses revealed that Guaymas geoducks 
are smaller and have a more inflated shell 
and deeper pallial sinus compared to those 
from San Felipe and Puerto Peñasco. Rocha- 
Olivares et al. (2010) also reported similar size 
differences between P. globosa from the upper 
and middle Gulf. 

Shape differences might also be explained by 
genetic differences, as observed in pectinids 
from the Pacific and Gulf sides of the Baja 
California Peninsula (Cruz & Ibarra, 1997; 
Cruz et al., 1998), or by phenotypic plasticity 
noted in other bivalve populations subjected 
to different environmental conditions (Palmer 
et al., 2004; Krapivka et al., 2007; Costa et 
al., 2008a; Elliot et al., 2008; Sokolowski et 
al., 2008; Márquez et al., 2010). In addition, 
morphological polymorphisms may result 
from gene-environment interactions (Rawson 
& Hilbish, 1991; Trapani, 2003). The variation 
in P. globosa we observed suggests either 
phenotypic plasticity or local adaptation due 
to gene-environment interaction. 

Shell shape differences between Bahía 
Magdalena and Gulf of California geoducks 
might be associated with temperate environ- 
mental conditions, in contrast to the northern 
Gulf, which is a warm-temperate region (Lluch- 
Belda et al., 2000; Calderón-Aguilera et al., 
2010). On the other hand, the lower shape 
differentiation between San Felipe and Puerto 
Peñasco might be caused by gene flow due 
to geographical proximity and circulation pat- 
terns, since the anticyclonic circulation during 
November-April in the northern Gulf (Lavín 
& Marinone, 2003) exactly overlaps with the 
reproductive period for P. globosa (Calderón- 
Aguilera et al., 2010). 

The notable shell shape differences of 
Panopea globosa in Guaymas may be due to 
drastically higher temperatures conditions dur- 
ing the summer months (Lluch-Cota & Arias, 
2000; Green et al., 2005); the most abundant 
beds on this site are located in shallower waters 
(8—12 m, Aragón-Noriega et al., 2007, Rocha- 
Olivares et al., 2010) compared with the other 
studied localities (Bahia Magdalena 10—20 m, 
González-Peláez Centro de Investigaciones 
Biológicas del Noroeste, México 2010, pers. 


comm.; San Felipe, Puerto Peñasco 10-25 
m, Rocha-Olivares et al., 2010). In addition, 
primary productivity is lower in Guaymas since 
there are essentially no currents most of the 
year (Lluch-Cota & Arias, 2000; Green et al., 
2005). 

The novel FB method of contours was found 
to have higher resolution than the traditional 
ILB method for distinguishing intraspecific phe- 
notypic differences among the three Panopea 
globosa groupings: (1) Bahía Magdalena; (2) 
San Felipe-Puerto Peñasco; and (3) Guaymas. 
Geometric morphometric analyses suggest the 
possibility of geographic and breeding isolation 
between FP. globosa of Bahia Magdalena in the 
Pacific and their conspecifics from the Gulf of 
California. Genetic studies are underway to un- 
derstand the evolutionary history of P generosa 
and P. globosa and the genetic relationships 
among populations of P. globosa. 


ACKNOWLEDGEMENTS 


Jesus Bautista and Carlos Pacheco of CIB- 
NOR provided technical help with sampling. 
Paul Valentich-Scott (Santa Barbara Museum 
of Natural History) provided valuable assis- 
tance with taxonomic identification. Delia Ro- 
jas-Posadas provided funding to obtain training 
in morphometrics. Ira Fogel, Laura Pérez and 
Miguel Cordoba of CIBNOR provided detailed 
editorial comments and constructive criticisms. 
We also thank Liza Ray and two anonymous 
reviewers for comments that greatly improved 
the manuscript. The research was funded by 
Consejo Nacional de Ciencia y Tecnologia 
(CONACYT, grant 106905). I.L.V. was the 
recipient of a SDN T doctoral fellowship 
(98412). 


LITERATURE CITED 


ALVAREZ, L., F. SUAREZ, R. MENDOZA & M. 
GONZALEZ, 2009, Bathymetry and active geo- 
logical structures in the upper Gulf of California. 
Boletín de la Sociedad Geológica Mexicana, 
61: 129-141. 

ANDERSON, L. & P. ROOPNARINE, 2005, Role 
of constraint and selection in the morphologic 
evolution of Caryocorbula (Mollusca: Corbul- 
lidae) from Caribean Neogene. Paleontologia 
Electronica, 8: 18 pp. 

ARAGON- NORIEGA, E., J. CHAVEZ, P. GRIB- 
BEN, E. ALCANTARA, A. MAEDA, E. ARAM- 
BULA, A. GARCÍA & R. MALDONADO, 2007, 
Morphometric relationships, gametogenic de- 
velopment and spawning of the geoduck clam 


MORPHOMETRICS IN PANOPEA 11 


Panopea globosa (BIVALVIA: HIATELLIDAE) in 
the central Gulf of California. Journal of Shellfish 
Research, 26: 423—431. 

BIZZARRO, J., 2008, A review of the physical 
and biological characteristics of the Bahia 
Magdalena lagoon complex (Baja California 
Sur, Mexico). Bulletin of the Southern California 
Academy of Sciences, 107:1-24. 

BOOKSTEIN, F., 1991, Morphometric tools for 
landmark data. Cambridge University Press, 
Cambridge, U.K., xvii + 435 pp. 

BRAY, N., 1988, Thermohaline circulation in the 
Gulf of California. Journal of Geophysical Re- 
search, 93: 4993-5020. 

BRIGGS, J., 1974, Marine zoogeography. 
McGraw-Hill, New York, x + 475 pp. 

BRUSCA, R., 1980, Common intertidal inver- 
tebrates of the Gulf of California, 274 ed. The 
University of Arizona Press, Tucson, Arizona, 
xx + 513 pp., [8] leaves of ols. 

CALDERON- AGUILERA, E., E. ARAGON- 
NORIEGA, W. -REYES C. PANIAGUA, A. 
ROMO & V. MORENO, 2010, Reproduction of 
the Cortez geoduck Panopea globosa (Bivalvia: 
Hiatellidae) and its relationship with temperature 
and ocean productivity. Journal of Shellfish 
Research, 29: 135—141. 

CARVAJAL, A., P. CONDE & E. ROLÁN, 2005, 
Decomposing shell form into size and shape 
by geometric morphometric methods in two 
sympatric ecotypes of Littorina saxatilis. Journal 
of Molluscan Studies, 71: 313—318. 

CARVAJAL, A., J. GUERRA, E. ROLAN & Á. E. 
ROLÁN, 2006, An example of the application 
of geometric morphometric tools to the mor- 
phological diagnosis of two sibiling species in 
Nassarius (Mollusca, Prosobranchia). Iberus, 
24: 1—88. 

COAN, E. V., P. VALENTICH-SCOTT & F. R. 
BERNARD, 2000, Marine bivalve mollusks from 
Arctic Alaska to Baja California. Santa Barbara 
Natural History Museum, Santa Barbara, Cali- 
fornia, viii + 764 pp. 

CONRAD, T. A., 1849, Fossils from northwestern 
America, Mollusca. Pp. 723-728, pls 17-21, in: 
J. Dana, ed., Geology. United States Exploring 
Expedition during the years 1838—1842. Under 
the command of Charles Wilkes, U.S.N., Vol. 
10, Appendix 1: xii + 10-756 pp.; Atlas, Geology, 
6 pp., 21 pls. 

COSTA, C., J. AGUZZI, P. MENESATTI, F. AN- 
TONUCCI, V. RIMATORI & M. MATTOCCIA, 
2008a, Shape analysis of different populations 
of clams in relation to their geographical struc- 
ture. Journal of Zoology, 276: 71—80. 

COSTA, C., P. MENESATT J. AGUZzZ!, S. 
D’ANDREA, F. ANTONUCCI, V. RIMATORI, 
F. PALLOTINO & M. MATTOCCIA, 2008b, 
External shape differences between sympatric 
populations of commercial clams Tapes decus- 
satus and T. philippinarum. Food Bioprocess 
Technol. Online: http://Awww.springerlink.com/ 
content/v37235277m7u8456/ 

CRUZ, P. &A. IBARRA, 1997, Larval growth and 
survival of two catarina scallop (Argopecten 
circularis, Sowerby, 1835) population and their 


reciprocal crosses. Journal of Experimental 
Marine Biology and Ecology, 212: 95—110. 

CRUZ, P., J. RAMIREZ, A. GUY & A. IBARRA, 
1998, Genetic differences between two popu- 
lations of catarina scallop (Argopecten ventri- 
cosus) for adaptations for growth and survival 
in a stressful environment. Aquaculture, 166: 
321-335. 

DALL, W., 1898, Contribution to the Tertiary 
fauna of Florida. Transactions of the Wagner 
Free Institute of Science of Philadelphia, 3(4): 
827-832. 

DEMERE, T. & N. SCOTT, 2006, Invertebrates 
of the Imperial Sea. Cap. 3, pp. 62—63, in: J. T. 
JEFERSON & L. Linpsay, eds., 2006, Fossil trea- 
sures of the Anza-Borrego Desert. Sunbelt Pub- 
lications, San Diego, California, xxi + 394 pp. 

DFO, 2000, Geoduck clam. Department of Fisher- 
ies and Oceans, Fisheries and Oceans, Canada. 
Science Stock Status Report C6-05: 3 pp. 

ELLIOT, J., K. HOLMES, R. CHAMBERS, K. 
LEON & P. WIMBERGER, 2008, Differences in 
morphology and habitat use among the native 
mussel Mytilus trossulus, the non native M. gal- 
loprovincialis, and their hybrids in Puget Sound, 
Washington. Marine Biology, 156: 39-53. 

FELDMAN, K., B. VADOPALAS, D. ARMSTRONG, 
C. FRIEDMAN, R. HILBORN, K. NAISH, J. 
ORENSANZ, J. VALERO, J. RUESINK, A. 
SUHRBIER, A. CHRISTY, D. CHENEY & J. DA- 
VIS, 2004, Comprehensive literature review and 
synopsis of issues relating to geoduck (Panopea 
abrupta) ecology and aquaculture production. 
Washington State Department of Natural Re- 
sources, Olympia, Washington, 140 pp. 

FITCH, J., 1953, Common marine bivalves of 
California. California Department of Fish and 
Game Marine Fisheries Branch, Sacramento, 
California, Fish Bulletin, 90: 102 pp., 1 pl. 

FYFE, A., 1984, The effect of conspecific associa- 
tion on growth and dispersion of the geoduck 
clam, Panope generosa. M.S. Thesis, Simon 
Frazer University, Vancouver, B.C., Canada. 

GOODWIN, L. & B. PEASE, 1989, Species pro- 
files: life histories and environmental require- 
ments of coastal fish and invertebrates (Pacific 
Northwest) — Pacific geoduck clam. U.S. Fish 
and Wildlife Service Biological Report No. 82 
(11.120). U.S. Army Corps of Engineers, TR 
EL-82-4: 23 pp. 

GOULD, A., 1850, Shells from the U.S. Exploring 
Expedition. Proceedings of the Boston Society 
of Natural History, 3: 214—218, 

GREEN, C., J. RUELAS & F. PAEZ, 2005, Mer- 
cury in surface sediments and benthic organ- 
isms from Guaymas Bay, east coast of the Gulf 
of California. Environmental Geochemistry and 
Health, 27: 321-329. 

GUMIEL, M., S. CATALA, F. NORIEAU, A. ROJAS 
DE ARIAS, A. GARCIA & J. DUJARDIN, 2003, 
Wing geometry in Triatoma infestatus (Klug) and 
T. melanosoma Martinez, Olmedo & Carcavallo 
(Hemiptera: Reduviidae). Systematic Entomol- 
ogy, 28: 173-179. 

HASTINGS, P., 2000, Biogeography of the tropi- 
cal eastern Pacific: distribution and phylogeny 


12 LEYVA-VALENCIA ET AL. 


of chaenopsid fishes. Zoological Journal of the 
Linnean Society, 128: 319—335. 

HELENES, J. & A. CARRENO, 1999 , Neogene 
sedimentary evolution of Baja California in 
relation to regional tectonics. Journal of South 
American Earth Sciences, 12: 589-605. 

INNES, D. & A. BATES, 1999, Morphological 
variation of Mytilus edulis and Mytilus trossulus 
in eastern Newfoundland. Marine Biology, 133: 
691-699. 

KEEN, M., 1971, Sea shells of tropical west 
America. Marine Mollusks from Baja California 
to Peru, 2nd ed. Stanford University Press, Stan- 
ford, California, xiv + 1064 pp., 22 pls. 

KESSLER, W., 2006, The circulation of the 
eastern tropical Pacific: a review. Progress in 
Oceanography, 69: 181-217. 

KONDO, Y., 1987, Burrowing depth of infaunal 
bivalves — observation of living species and its 
relation to shell morphology. Transactions of the 
Proceedings of the Palaeontological Society of 
Japan, n.s., 148: 306-323. 

KRAPIVKA, S., J. TORO, A. ALCAPÁN, M. AS- 
TORGA, P. PRESA, M. PÉREZ & R. Guiñez, 
2007, Shell- -shape variation along the latitudi- 
nal range of the Chilean blue mussel Mytilus 
chilensis (Hupe 1854). Aquaculture Research, 
38: 1770-1777. 

KWON; J., J. PARK, Y. LEE, J. PARK, Y. HONG 
& Y. CHANG, 1999, Morphological variation and 
genetic relationship among populations of short 
necked clam Ruditapes philippinarum collected 
from different habitats. Journal of Fisheries Sci- 
ence and Technology, 2: 98—104. 

LAVIN, M. & S. MARINONE, 2003, An overview 
of the physical oceanography of the Gulf of 
California. Pp. 173—204, in: O. VELASCO-FUENTES, 
J. SHEINBAUM & J. Ochoa, eds., Nonlinear pro- 
cesses in geophysical fluid dynamics. Kluwer 
Academic Publishers, Dordrecht, Netherlands, 
XP O7Gipp. 

LAWING, A. & P. POLLY, 2009, Geometric mor- 
phometrics: recent applications to the study of 
ae and development. Journal of Zoology, 

80: 1-7. 

LOCKWOOD, R., 2004, The K/T event and infau- 
nality: Morphological and ecological patterns of 
extinction and recovery in veneroid bivalves. 
Paleobiology, 30: 507-521. 

LLUCH-COTA, S. & J. P. ARIAS, 2000, Sobre 
la importancia de considerar la existencia de 
centros de actividad biológica para la regional- 
izacion del océano: el caso del Golfo de Cali- 
fornia. Pp. 255-263, in: D. LLucH, J. ELOURDUY, 
S. LLucH & G. Ponce, eds., BAC. Centro de 
Investigaciones Biológicas del Noroeste, La 
Paz, B.C.S., Mexico, 362 pp. 

LLUCH- BELDA, D., M. HERNANDEZ, R. 
SALDIERNA & R. GUERRERO, 2000, Vari- 
abilidad de la temperatura superficial del mar 
en Bahía Magdalena, B.C.S. Oceánides, 15: 
1—23. 

LLUCH-COTA, S., E. ARAGÓN-NORIEGA, F. 
ARREGUÍN, D. AURIOLES, J. BAUTISTA, R. 
BRUSCA, R. CERVANTES, A. CORTES, P. DEL 
MONTE, A. ESQUIVEL, G FERNANDEZ, M. 


HENDRICKX, S. HERNANDEZ, H. HERRERA, 
M. KAHRU, M. LAVIN, D. LLUCH- BELDA, D. 
LLUCH- COTA, di LOPEZ, S. MARINONE, M. 
NEVAREZ, S. ORTEGA, E. PALACIOS, A. 
PARES, G. PONCE, M. RAMIREZ, Cc. SALINAS, 
R. SCHWARTZLOSE &A. SIERRA, 2007, The 
Gulf of California: review of ecosystem status 
and sustainability challenges. Progress in 
Oceanography, 73: 1—26. 

MAHALANOBIS, P., 1936, On generalized dis- 
tance in statistics. Proceedings of the National 
Institute of Sciences of India, 2: 49-55. 

MARKO, P. & J. JACKSON, 2001, Patterns of 
morphological diversity among and within 
arcid bivalve species pairs separated by the 
Isthmus of Panama. Journal of Paleontology, 
75: 590-606. 

MARQUEZ, F., J. ROBLEDO, G. ESCATI & S. 
VAN DER MOLEN, 2010, Use of different geo- 
metric morphometrics tools for the discrimina- 
tion of phenotypic stocks of the striped clam 
Ameghinomya antiqua (Veneridae) in north 
Patagonia, Argentina. Fisheries Research, 101: 
127—131. 

MATIAS, A., J. DE LA RIVA, M. TORREZ & J. 
DUJARDIN, 2001, Rhodnius robustus in Bo- 
livia, identified by its wings. Memories Institute 
Oswaldo Cruz, 96: 947—950. 

MOORE, E., 1968, Fossil mollusks of San Diego 
County. San Diego Society of Natural History, 
Occasional Paper, 15: 76 pp. 

MOORE, T., 2001, Geoduck. Pp. 449—450, in: W. 
Leet, CH. Dewees, R. KLINGBEIL & E. LARSON, eds., 
California’s living marine resources: a status 
report. California Department of Fish and Game. 
Resources Agency. http://portal.nceas.ucsb. 
edu/working_group/science-frameworks-for- 
ebm/background-information/CA Living Marine 
Resources statusreport.pdf 

MORRIS, R., D. ABBOT & E. HADEDLE, 1980, 
Intertidal Invertebrates of California. Stanford 
University Press, Stanford, California, 388 pp. 

PALMER, M., G. PONS & M. LINDE, 2004, Dis- 
criminating between geographical groups of a 
Mediterranean commercial clam [Chameleagal- 
lina (L.): Veneridae] by shape analysis. Fisheries 
Research, 67: 93-98. 

RAWSON, P. & T. HILBISH, 1991, Genotype- 
environment interaction for juvenile growth 
in the hard clam Mercenaria mercenaria (L.). 
Evolution, 5: 1924—1935. 

ROCHA-OLIVARES, A., L. CALDERON-AGUIL- 
ERA, E. ARAGON- NORIEGA, N. SAAVEDRA & 
V. MORENO, 2010, Genetic and morphological 
variation of northeast Pacific Panopea clams: 
evolutionary implications. Journal of Shellfish 
Research, 29: 327—335. 

ROHLF, F. J., 1990, Morphometrics. Annual Review 
of Ecology and Systematics, 21: 299—316. 

ROHLF, F. J., 2008, tos Dig2. http://life.bio.sunysb. 
edu/morph/index.html 

ROHLF, F. J., 2009, tps Util Version 1.44. http:// 
life.bio.sunysb.edu/morph/soft-utility.html 

ROOPNARINE, P. & G. VERMEIJ, 2000, One 
species becomes two: The case of Chione 
cancellata, the resurrected C. elevata, and a 


MORPHOMETRICS IN PANOPEA 13 


phylogenetic analysis of Chione. Journal of 
Molluscan Studies, 66: 517-534. 

ROOPNARINE, P., J. SIGNORELLI & C. LAUM- 
ER, 2008, Systematic, Biogeographic and 
microhabitat-based morphometric variation of 
the bivalve Anomalocardia squamosa (Bivalvia: 
Veneridae: Chioninae) in Thailand. The Raffles 
Bulletin of Zoology, Suppl. 18: 95—102. 

SAVAZZI, E., 1987, Geometric and functional 
constraints on bivalve shell morphology. Leta- 
hia, 20: 293-306. 

SAGARPA, 2007, Secretaria de Ganaderia, Agri- 
cultura, Desarrollo Rural, Pesca y Alimentacion, 
Mexico. Informe de Labores. www.sagarpa.gob. 
mx/infohome/pdf/labores_sagarpa.pdf 

SHEETS, D., C COVINO, J. PANASIEWICZ & 
S. MORRIS, 2006a, Comparison of geometric 
morphometric outline methods in the discrimina- 
tion of age-related differences in feather shape. 
Frontiers in Zoology, 3: 15. 

SHEETS, D., M. ZELDITCH & D. SWIDERSKI, 
2006b, IMP (Integrated Morphometrics Pack- 
age): http://www.canisius.edu/~sheets/morph- 
soft.html 

SHIPUNOV, A. & R. BATEMAN, 2005, Geomet- 
ric morphometrics as a tool for understanding 
Dactylorhiza (Orchidaceae) diversity in Euro- 
pean Russsia. Biological Journal of the Linnean 
Society, 85: 1-12. 

SOKOLOWSKI, A., K. PAWLIKOWSKI, M. 
WOLOWICZ, P. GARCIA & J. NAMIESNIK, 
2008, Shell deformations in the Baltic clam 
Macoma balthica from southern Baltic Sea (the 
Gulf of Gdansk): hypotheses on environmental 
effects. Ambio, 37: 93—100. 

STANLEY, S., 1970, Relation of shell form to life 
habits of the Bivalvia (Mollusca). Geological 


Society of America Memoir, 125: xiii + 296 pp., 
40 pls. 

STANLEY, S., 1979, Macroevolution: pattern 
and process. Johns Hopkins University Press, 
Baltimore, Maryland, xi + 332 pp. 

TRAPANI, J., 2003, Morphological variability in 
Cichlasoma minckleyi, the Cuatro Cienegas 
cichlid. Journal of Fish Biology, 62: 276—298. 

UBUKATA, T., 2003, A morphometric study on 
morphological plasticity of shell form in crevice- 
dwelling Pterioida (Bivalvia). Biological Journal 
of the Linnean Society, 79: 285—297. 

VADOPALAS, B., T. PIETSCH & C. FRIEDMAN, 
2010, The proper name for the geoduck: resur- 
rection of Panopea generosa Gould, 1850, from 
the synonymy of Panopea abrupta (Conrad, 
1849) (BIVALVIA: MYOIDA: HIATELLIDAE). 
Malacologia, 52: 169-173. 

WEYMOUTH, F., 1921, The edible clams, mussels 
and scallops of California. State of California 
Fish and Game Commission, Fish Bulletin, 4: 
74 pp., 19 pls. 

YONGE, C., 1971, On Functional morphology 
and adaptative radiation in the bivalve super- 
family Saxiscavacea (Hiatella (= Saxicava), 
Saxicavella, Panomya, Panope, Cyrtodaria). 
Malacologia, 11: 1—44. 

ZAYTSEV, O., R. CERVANTES, O. MONTANTE & 
A. GALLEGOS, 2003, Coastal upwelling activity 
on the Pacific shelf of the Baja California Penin- 
sula. Journal of Oceanography, 59: 489—502. 

ZELDITCH, M., D. SWIDERSKI, D. SHEETS & 
W. FINK, 2004, Geometric morphometrics for 
biologists: a primer. Elsevier/Academic Press, 
San Diego, California, x + 443 pp. 


Revised ms. accepted 12 July 2011 


